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Abstract 

The Cell Network Model is a fracture model recently in- 
troduced that resembles the microscopical structure and 
drying process of the parenchymatous tissue of the Bam- 
boo Guadua angustifolia. The model exhibits a power-law 
distribution of avalanche sizes, with exponent —3.0 when 
the breaking thresholds are randomly distributed with uni- 
form probability density. Hereby we show that the same 
exponent also holds when the breaking thresholds obey a 
broad set of Weibull distributions, and that the humidity 
decrements between successive avalanches (the equivalent 
to waiting times for this model) follow in all cases an ex- 
ponential distribution. Moreover, the fraction of remaining 
junctures shows an exponential decay in time. In addition, 
introducing partial breakings and cumulative damages in- 
duces a crossover behavior between two power-laws in the 
avalanche size histograms. This results support the idea 
that the Cell Network Model may be in the same univer- 
sality class as the Random Fuse Model. 

Statistical models of fracture Finite Element Method 
Computational mechanics of solids. 

PACS 02.50.-r,05.90.+m,46.50.+a,62.20.F-, 62.20.M- 

1 Introduction 

The Cell Network Model of Fracture (CNMF [I]), is a two 
dimensional statistical model of fracture [2] inspired in the 
stress field caused by drying of the bamboo Guadua angus- 
tifolia [2111]. At the parenchymatous tissue level, bamboos 
shrink during drying, causing the detaching of neighboring 
cells and the appearance of fractures. The CNMF mod- 
els this tissue as an hexagonal array of cell elements (each 
of them made of six beams), fixed by angular springs and 
joined by brittle springs called the junctures. (Figure: [l]). 
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Shrinking forces -due to drying- acting along the elements 
distort the structure and cause breaking avalanches of the 
junctures among cells. 




Figure 1: Cell Network Model of Fracture (CNMF). Upper 
left, The plane frame element spanning between the nodes 
i and j, oriented by an angle 9. Each node has two trans- 
lational and one rotational degrees of freedom. Lower left. 
Two contiguous Cells. Right Structure of the CNMF. The 
hexagons represent the cells and the junctures are arranged 
into triangles. Hashing denote fixed boundary conditions. 
From pQ Not at scale. 

Interestingly enough, when an homogeneous distribution 
of breaking thresholds of junctures and fixed Young mod- 
ulii are used, the histogram of avalanche sizes of the CNMF 
shows power law behavior with an exponent of —2.93(8) 
(Q]). This is by all means equal to avalanche size distribu- 
tion of the random fuse model, which shows a power law 
with exponent of —3 [3[6]. The analytical solution of both 
models has been elusive. 

One dimensional fiber bundle systems have provided 
models for the thoroughly study of critical phenomena. 
Universality, the effect of damage and the critical expo- 



nents have been found analytically (see [7[[S] and references 
therein.). The logical extension of those models to 2D, the 
random fuse model, has allowed to investigate the fracture 
properties of biological materials, as in the case of brittle 
nacre, ([!]), describing the toughness of the material by 
its microscopical architecture. Beam models similar to the 
CNMF has also been used to model fracture in concrete 

In the present paper, to characterize the path to global 
failure, we study both the distribution of humidity decre- 
ment among between consecutive avalanches (the analogue 
of waiting times for this model), as well as the fraction of 
intact fibers as function of the humidity decrement. This 
last quantity behaves as an order parameter for the system 
(see (Sec: ^jj). Moreover, the universality of the CNMF 
is explored numerically by characterizing the histogram of 
avalanche sizes for two cases: homogeneous distributions of 
the juncture breaking thresholds with different widths and 
Weibull distributions of different shapes, (Sec: [3]). Fur- 
thermore, in section (Sec: Q a damage function is intro- 
duced as follows: Each time a stress threshold is reached, 
the stiffness is reduced by a constant factor, until the fiber 
completely breaks. The main results and comments are 
summarized in Sec (Sec: [5]). 

2 Model 

The CNMF is a 2D statistical model of fracture that re- 
sembles the parenchymatous tissue of the bamboo Guadua 
augustifolia. It is composed by two kinds of structures: 
cell walls and junctures among cells. Six cell walls arrange 
themselves to form hexagons, thanks to a angular springs 
associated with the rotational degree of freedom. The cells 
are arranged like a honeycomb. The junctures are arranged 
in sets of three at the common corners of the cells, mod- 
eling the silica deposits that glue cells together. (Figure: 
[TJ Each kind has a given fixed Young modulus for all its 
elements. Junctures are allowed to break, as a result of the 
brittle behavior of the silica, while cell walls are not. 

As boundary conditions, all the border nodes are set 
fixed. The deformation of the system comes from shrink- 
ing forces acting on every cell wall and proportional to a 
global humidity loss parameter Ah. By means of a Finite 
Element Method, the resulting forces and deformations of 
all the elements are calculated. 

The evolution of the system has three stages. Linear 
Elastic Shrinking: local shrinking forces due to humidity 
losses are applied to the cell walls. The differences between 
the local strain at the juncture and their individual thresh- 
olds are calculated. This step continues until at least one 
fiber would suffer an strain surpassing its threshold. Dry- 
ing induced Breaking: By means of a zero finding algo- 
rithm, the exact humidity loss causing the first breaking 
is found. The broken element is removed from the struc- 



ture, changing the stiffness matrix in accordance. Nonlin- 
ear avalanche: The breaking of an element calls for a force 
redistribution over the whole structure. This redistribution 
may cause an avalanche of breaks. When the avalanche 
ends, the procedure re-stars from the first stage. 

2.1 Distributions of humidity decrement 
between consecutive avalanches 

The distribution of humidity decrements between succes- 
sive avalanches provides a description of the path to the 
global failure of the system. This is the analog to the wait- 
ing times between avalanches of other models of fracture. 
Fig. (Figure: [2]) shows the normalized histograms of hu- 
midity decrements for several values of the maximum hu- 
midity change allowed Ah. When the breaking process is 
driven until the end, the histograms can be fitted to an ex- 
ponential, with a fitted humidity constant of 13.5(2), which 
is and indication of lack of correlation between successive 
avalanches. 




Figure 2: Humidity decrement between consecutive 
avalanches histogram. The parameter that defines a curve 
is the maximum humidity change (Ah, proportional to Cell 
Wall strain, shown). The classes (horizontal axis) are the 
(normalized by the mean waiting humidity). 



2.2 Fraction of remaining junctures 

Let us consider the one-dimensional global load sharing 
fiber bundle model |S]. Given U(a) as the fraction of 
remaining fibers at a given stress and o~ c the critical 
stress causing the global breaking of the system. Thus, 
U*(a) — U*(<j c ) behaves as an order parameter. For the 
CNMF with homogeneous breaking thresholds at the junc- 
ture elements we obtain an exponential relaxation in the 
number of remaining junctures ((Figure: [3])). 
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Figure 3: Number of intact fibers as function of the hu- 
midity change for several homogeneous distribution of the 
breaking threshold with different widths (semi log). Inset, 
linear axis. 



This exponential behavior is independent of the system 
size. In (Figure: [2]), the fraction of the population of intact 
fibers as function of humidity is shown for different sizes of 
the system. The horizontal axis was rescaled by means of 
a linear fit on the semi-log data. All histograms show an 
exponential decay. 



teristics. For the CNMF we studied the distribution of 
avalanche sizes when the thresholds are generated either 
from flat distributions of several widths or from Weibull 
distributions with several characteristic parameters. 

Flat distributions of the breaking thresholds, all centered 
at 0.35EA but with different widths (spanning on two or- 
ders of magnitude), show the same power-law distribution 
of avalanche sizes, with slopes around —3 (Fig. (Figure: [5|). 
The data shows that narrower distributions show larger 
fluctuations around this power law than wider ones. 
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slope = -3.0(1) 
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Figure 4: Scaled fraction of intact fibers as function of the 
humidity change (in units of max humidity difference) for 
different system sizes. 



3 Universality 

Fiber bundle models are universal in the sense that the 
breaking of the elements follows power law distribution 
of avalanche sizes irrespective of several system charac- 



Figure 5: Histogram of avalanche sizes for several widths 
of the homogeneous threshold distribution centered at 
0.35EA 



The probability density function for the Weibull distri- 
bution function [11] is given by: 



f(x;X,k) 



(!) 



k-i 



x>0 
x < 



(1) 



where k and A are free parameters. It is commonly used 
to describe the breaking thresholds of fibers by fixing A=l 
and changing fc, which is the main parameter controlling 
the distribution shape. When we use this distribution for 
the breaking thresholds at the junctures, the histogram of 
avalanche sizes shows also a power law behavior, with an 
exponent close to -2.9 for small values of k. For larger 
values the exponential cutoff is more pronounced. (Figure: 

©• 

The humidity decrements between consecutive 
avalanches (that is, the waiting times) distributes 
like an exponential, also when the thresholds follow a 
Weibull distribution ((Figure: [7|).The characteristic time 
for k=5 is 16.1(8), very close to the one we gathered for 
flat distributions of the braking thresholds. 
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Figure 6: Histogram of avalanche sizes for several Weibull 
distributions of the breaking thresholds, with values of the 
shape parameters k between 1 and 10. The straight line 
corresponds to an slope of —2.9. 

4 Damage 

In order to introduce a degradation for the juncture ele- 
ments we reduce the Young modulus of the juncture ele- 
ments by a damage factor < a < 1 each time the juncture 
fails. When the element has suffered a maximum number 
of failures k max , it is assumed to be broken and is removed 
from the structure. 

When a small degradation is introduced (a— 0.1), the 
power law distribution of avalanche sizes seems to exhibit a 
crossover from an exponent —3 to an exponent —2 at sizes 
around 8 (Fig. (Figure: [8|). This may indicate that the 
remaining elasticity helps to sustain the structure. How- 
ever, more statistics and larger system sizes are required 
to clarify this point. Even smaller values of a (not shown) 
cause the maximum humidity loss (and therefore the forces 
on the elements) to be much larger, creating numerical in- 
stabilities that end into poor statistics. 



5 Conclusions 

The numerical evidence of this work indicates that the 
avalanche sizes for the Cell Network Model of Fracture dis- 
tribute as a power law with exponent —3.0, for any broad 
distribution of the braking thresholds, either if they are flat 
or Weibull distributed. 

The distribution of waiting times show an exponential 
decay for all system sizes evaluated and all tested disorder 
distributions of breaking thresholds. Even the characteris- 
tic times are similar for all of them. In our opinion, this is 
related to the fact that the most common failure mode for 
the system is the softening of the sample by displacements 
that would violate the boundary conditions. 



Figure 7: Histogram of humidity increments between suc- 
cessive avalanches for several Weibull distributions of the 
breaking thresholds, with k between 1 and 10. The blue 
line corresponds to a fit of the series of shape k = 10, with 
slope —2.1. The green line to that for k = 3, with slope 
-2.4 
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